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Abstract 

^— ( ' We study the renormalisation group running of the cosmological and the Newton constant, where 

the renormalisation scale is given by the inverse of the radius of the cosmological event horizon. In this 
CO ' framework, we discuss the future evolution of the universe, where we find stable de Sitter solutions, 

^ ' but also "big crunch"-like and "big rip"-like events, depending on the choice of the parameters in the 

00 model. 

. 

o ■ 

O ; 1 introduction 

o ■ 

■ The recently observed [1] accelerated expansion of the universe may have its reason in the existence of dark 
energy (DE), an energy form with negative pressure, which is so far not understood. The cosmological 

^ ■ constant (CC) is the theoretically simplest candidate for DE, because it occurs as a classical parameter 

in Einstein's equations, and further it has an origin as vacuum energy in quantum field theory (QFT). 

■ On the other hand, it is difficult to explain the tiny value of the CC and the actual coincidence of the 
. !^ ' energy densities of the CC and non-relativistic matter [2] . 

, There are lots of models, which describe DE as a dynamical quantity, e.g. by using scalar fields. Another 

' possibility is the modification of the theory of gravity by introducing extra terms in the equations of the 

cosmological evolution, or extending our space-time by additional space-time dimensions. However, in 
most of these models the accelerated cosmic expansion is due to new and unknown physics, which often 
means a high amount of arbitrariness and limited predictability. 

In this work, we investigate the CC in the sense that it emerges anyway on a formal level in QFT. There, 
the zero-point or vacuum energy of a quantised field has the same equation of state as the CC occurring 
in Einstein's equations. Unfortunately, it is unknown how to calculate its value in a unique way, because 
it can be written in the form of a quartically divergent momentum integral like J d^p ■ p. The naive 
assumption of an ultra-violet (UV) cutoff in this integral at some known energy scale usually leads to an 
unobserved high value of the CC, which is also called the old CC problem. However, the procedure of 
renormalisation of (coupling) constants in QFTs can handle infinities, thereby leading to a dependence 
of the renormalised constants on some energy scale /i. In many cases, this renormalisation scale can be 
identified with an external momentum, or at least with some characteristic scale (e.g. the temperature) 
of the environment. Studying QFT on curved space-time [3] leads to infinities in the vacuum expectation 
values (VEV) of the energy-momentum tensors of the fields. This can be treated by renormalisation to 
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yield a scale-dependent or running CC and a running Newton constant (NC). The absolute values are 
still not calculable, but the change with respect to the renormalisation scale can be described by the 
renormalisation group equations (RGE). Unlike the running coupling constants in the standard model 
of particles, here, the physical meaning of the scale is not given by the theory. This becomes obvious 
since, in the language of Feynman graphs, the formally infinite value of the CC corresponds to a closed 
loop without external legs and hence no distinct energy scale. Connecting this renormalisation scale with 
physics thus requires an additional theoretical input. Usual choices in the literature are the Hubble scale, 
the square root of the Ricci scalar, and combinations of similar quantities. The RG running of the CC 
and the NC has been studied in several different frameworks and models, some recent results can be 
found in Rcfs. [4-8]. 

In our investigation of the RG running of the CC and the NC, we choose the inverse of the radius of 
the cosmological event horizon as the renormalisation scale. Such a horizon usually exists in accelerating 
universes like ours. In addition, the possible relevance of the horizon scale for DE has often been pointed 
out [9-11]. In section 13 we derive the RGEs and their dependence on the renormalisation scale and 
the parameters, mainly the field masses. Since the event horizon in an evolving universe is not constant 
in time, the CC and the NC are also time-dependent, implying that the usual covariant conservation 
equations for the energy momentum tensor have to be modified (section EJ. In section 0] we discuss the 
properties of the new evolution equation for the cosmic scale factor and derive some conditions on the 
existence and the stability of the final states of the universe. The cosmic fate is the main point of the 
discussion, since the characteristic behaviour in the far future depends crucially on the parameters in 
the RGEs. In sectional we illustrate the possible final states of the universe by showing some numerical 
solutions and their dependence on the parameters. Finally, section El contains our conclusions and some 
open points of this setup. 



2 Renormalisation group equations 



To formulate the RGEs for the CC and the NC G, we consider free quantum fields on a curved space-time, 
namely a Friedmann-Robertson- Walker (FRW) universe with a positive CC A. For one fermionic and 
one bosonic degree of freedom with masses mp and tob, respectively, the 1-loop effective action can be 
written in the form [3] 



ScS = I d^a;^ 



+ C, (1) 



where D = ^7Euicr + lim„^4(n — 4)~^ is a divergent term, which does not depend on the renormalisation 
scale /i. Furthermore, ^ is a coupling constant^, and the variable C represents all further terms in 
the effective action, that are neither proportional to the Ricci scalar Ric nor to the vacuum energy 
density A A/(87rG). 

The relevant /3-functions in the MS-scheme for the vacuum energy density A and the NC G are obtained 
by the requirement that the effective action S^s must not depend on the renormalisation scale /z, 

dS'off „ 

0. 

d/i 

Because of this condition, A and G have to be treated as //-dependent functions in Eq. which have 
to obey the RGEs given by 



dA mt — mi, d / 1 \ 1 r/, i 



d/i 327r2 dfi \G J tt 



Note, that the divergent term D has dropped out, leaving over just the masses mp/B and ^. Assuming 
constant masses, the RGEs can be integrated. Hence, the equation for the vacuum energy density reads 

A(/.) ^ Ao(l - gi In^), Ao:=A(/io), (2) 



-"^In the action S = J d^x^—g ^ ^j'g^g'^ -I- ^4':a4"°' ~ fi™^ -I- 5R-ic](/'^j of a scalar field on a curved space-time, the 
constant ^ occurs in the coupling term ^ ■ Ric • (p'^ between the scalar field and the Ricci scalar Ric. 
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where the sign of the parameter 

'^^ 32^0 ~ 

depends on whether bosons or fermions dominate. In this context, a real scalar field counts as one bosonic 
degree of freedom, and a Dirac field as four fermionic ones. The generalisation to more than one quantum 
field in the RGE, can be achieved by summing over the fourth powers of their masses. For the NC G we 
obtain the RGE in the integrated form 

= , ^° , , Go:=G(/io). (4) 
1 - 52 In - 

Again, we omit the generalisation to more fields, that follows from summing over the squared masses of 
the fields. For one bosonic and one fermionic degree of freedom the mass parameter q2 is given by 

G 

12 [{i - ^)ml - j^ml] . (5) 

Finally, we remark, that Eq. (0) for the running vacuum energy density A(/i) was derived in a renormali- 
sation scheme, which is usually associated with the high energy regime. Unfortunately, the corresponding 
covariantly derived equations for the low energy sector are not known yet [12]. Therefore, we prefer to 
work with the above RGEs, which were derived in a covariant way, and study the consequences and the 
constraints on the mass parameters qi and (72 • 

Next, we choose the renormalisation scale 11 to be the inverse of the radius R of the cosmological event 
horizon. In the FRW universe the (radial) horizon radius R at the cosmic time t is given by 

r°° dt' 

H-' = R{t):^a{t)- — -, (6) 
Jt ait') 

where a(t) is the cosmic scale factor, corresponding to the line element^ 

ds^ = df - a^{t) ( + r^dn' 
\ 1 — Kr^ 

For universes, which end within finite time, the upper limit of the integral in Eq. should be replaced 
by the time when the universe ends. 

The choice of the scale fi = R^^ can be motivated by the thermodynamical properties of the cosmological 
event horizon. This horizon emits radiation, whose temperature is given by the Gibbons-Hawking [13] 
temperature Tgh = {27ri?)~^, that is proportional to our renormalisation scale /i. For a comoving observer 
in a de Sitter universe the only cosmological energy scale is given by this temperature. We have to admit 
that this no prove of the rightness of this choice. On the other hand, the investigation of the cosmological 
evolution with this specific renormalisation scale is the main point of this work, and the resulting solutions 
are quite interesting. From Eq. JBJ we see that in an evolving universe the event horizon radius R and 
thus the scale /j, are usually not constant in time. Therefore, the vacuum energy density A and the NC G 
will be time-dependent, too. This requires a generalisation of the covariant conservation conditions, and 
complicates the solutions of Friedmann's equations. 



3 Evolution equation for the scale factor 

In this section, we derive the evolution equation for the cosmic scale factor a{t) in the framework of 
the spatially isotropic and homogeneous FRW universe with a time-dependent CC and NC. On this 
background, radiation and pressureless matter (dust) can both be described by a perfect fluid with the 
energy density p and the pressure p = cup, where the constant cu characterises the equation of state"^. The 
corresponding energy-momentum tensor for these energy forms reads 

T"'' = (p + p)u"w'3-ra"^, 

■^The constant k fixes the spatial curvature of the universe. 
''Dust: LO = 0; radiation: u) = |. 
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with being the four-velocity vector field of the fluid. With our choice of the renormalisation scale, G 
and A depend only on the cosmic time t. From Einstein's equations 



conservation equations'* 



;/3 



and from the contracted Bianchi identities G"^g = for the Einstein tensor G"^, we obtain the generalised 



= [GA5"'' + Gr"^];/3 = G(A + p) + G(A + p + 3fp(l + c^)). 

Note that the simple scaling rule p (x a~^^*+'^^ for the matter content is not valid anymore, because it is 
now possible to transfer energy between the matter and the vacuum, in addition to G 7^ 0. Therefore, 
we have to combine the Friedmann equations for the Hubble scale H := ^ and the acceleration ^, 

+ ^ - yGW(A(0 + pW), (7) 

- = ^G{t){A{t)~Qp{t)), Q:=i(l + 3c.), (8) 
ad 2 

to eliminate the matter energy density p. The left-hand side of the result is abbreviated by F{t): 



m := - 



2 



k 



= |^GA.(l + g). (9) 



Now the RGEs for A and G from Eqs. Q and Q are inserted in Eq. © to yield, with the specific choice 
of our renormalisation scale p = 1/R, the main equation of this work: 



Mo R{t) 

— = exp 



p{t) Ro 

In this equation the constant Kq is defined as 



KoF{t) - 1 



qi - q2KoF{t) 



(10) 



3 H^^ 



87rGoAo(Q + l) f^AoCQ + l)' ^^^^ 

where 17 ao = 87rGoAo/(3-ffo) the relative vacuum energy density and Hq the Hubble scale at the 
time t = to- 

Note that for a dominant matter energy density p^ A and flat spatial curvature (fc = 0), the acceleration 
quantity g := ^ is given by the negative value of the (new) equation of state parameter Q = (1 + 3cj)/2, 
which we have introduced in Eq. 0. 



4 Discussion 

This section is devoted to the discussion of the properties of our main equation (|10|l . thereby placing 
special interest in the late-time behaviour of the scale factor a{t). Solving Eq. (|10|l for F leads to 



1 + gi In 

KoF^-——^, (12) 
1 + 92 In j| 



which is plotted as a function of the horizon radius R in Figs. 01 El El and for four cases, depending 
on the signs of qi and q2. For reasonable fleld masses m, the magnitude of the mass parameter qi is 
always greater than that of q2. In the definition of the parameter q2 the field masses m are divided by 
today's Planck mass Afpianck.o = 1 / V Go , therefore, |(72| is a tiny quantity today. This implies a very weak 
running of the NC G, which agrees with the strong bounds on the time- variation of G [14]. Additionally, 
this has the advantage, that today we are far away from the Landau pole of G{p), where the function F 



■'We do not assume T' 



,al3 
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Figure 1: Left: The real-valued branches of Lambert's VK-function Wu{z), u — 0, —1. Right: The lower 
bound -Romin of the initial event horizon radius Ro for negative gi. For Ro < i?omin a final de Sit- 
ter state does not exist. 



diverges. Very high values of |i^| could render our calculations invalid, since we do not implement higher 
powers of the curvature scalar in the gravitational action. These contributions are probably relevant in 
regimes with large Whether they are able to prevent the singular behaviour in the scale factor, that 
occurs in the numerical solutions for some cases, requires further investigation and is not treated in this 
paper. 

Much less problematic are final states of the universe, that are de Sitter-like. In this case, the scale 
factor acquires the form a{t) oc exp{Hct) for large t, where He denotes the constant Hubble scale, and the 
radius of the event horizon is given by i?o = l/^fc- Plugging this asymptotic form for a{t) in Eq. (|12|l . 
one arrives at 

Ko{l + Q) _2 l + gilnx 

/to 1 + 92 in a; 

where the variables x :— Rc/Rq and Ko{l + Q)/Rq > have been introduced. Neglecting the 

running of the NC, we set q2 = and therefore we have to solve qax^^ = 1 + qi Ina; for x. The results 
are given by 

""^To^ ^ q,.Wu{^e^^y ^^^^ 

involving Lambert's VF-function Wu{z), which is the inverse function of z = xe^ . The index u denotes the 
different branches of this function. Only for u = 0, — 1 it takes on real values for real arguments z > — e^^. 
Additionally, W-i{z) is not real- valued for 2 > 0, as can be seen in Fig.Q] where both branches are plotted. 
From these properties of the M^-function we get for negative qi the constraint qs < — ^exp(— ^ — 1), 
which implies an lower bound for the initial value Rq of the horizon radius: 



2 / 2 

Ra > itomin := \ 7Z — TH^^P ^ 

y qiilAoH^ \qi 

If Rq is smaller than this minimal value, then Eq. (|13|l has no positive solutions and a final de Sitter 
state does not exist. For Rq = i?omin there is exactly one solution x — exp(— ^ — i), for higher values Rq 
there are two solutions. In the case of a positive value of the parameter gi, the initial value Rq must be 
smaller than 1/ ^HqHao- Otherwise the final horizon radius Re is smaller than the initial one, a; < 1. 
Both cases are plotted in Fig. [3 

Since we have found several de Sitter solutions, we have to study the stability of these final states. 
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Figure 2: The ratio x between the final event horizon radius Rc and the initial radius Ro as a function of Ro, 
Eq. I13II . Here qi is set to ±2 and q2 — 0. All stable de Sitter final states lie on the curve q2 — —2, u — 0. 

Therefore, we write KqF as a function of KqF = Ko{H + {Q + 

KqF = qi - — exp 

where we used R = RH — 1. In the final de Sitter state we have KqF = and R = Rc = l/i?e = const. 
Near this point we can neglect H in the function F and replace H by ^ Kp^g+ iy ■ 
it is required that 

d(A'oi^) 



For a stable solution 



d{KoF) 



< 



in the final point, where KqF = Kq{Q + ^)F[^ 
stability condition 



q^x ^. With and x from above, this yields the 



< -1, 



implying that there are no stable de Sitter solutions for positive values of qi , because the VF-function is 
positive. For negative qi we get the condition 

W^C-^e'/"^) > -1, 

which means that only the solution with u = is stable. This renders the final event horizon radius Rc = 
Rqx unique. 

Finally, we take a closer look at the ratio Rc/Rq as a function of the mass parameter qi. For initial 
values i?o < l/(^foV^Ao)j which means ga > 1, there is a certain range of values of qi where no solutions 
for Rc exist. This range is again given by the requirement that the argument of the M^-function must be 
greater than or equal to — e~^, leading to the conditions 



91 < 



or qi > 



(14) 



In Fig. 131 the exclusion range for qi is obvious for q^ > 1. In the case that qi lies above this range, the 
unstable solution for Rc is reached first during the future cosmic evolution. For qi below this range, the 
stable solution is nearer to the initial value Rq than the unstable one, however, both solutions for Rc lie 
below Rq. Initial values Rq > l/(7Jo\/f^Ao) (i-c 93 < 1) lead to stable final states with Rc > Rq for all 
negative values of qi . 
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Figure 3: This figure shows the ratio Rc/Ro between the final and the initial event horizon radius as a function 



of qi. We plotted four cases for different values of gs = l/{HgRQQAo), where the choice gs = 0, 8; 0, 9 
corresponds to an initial radius Ro > l/(-ffo\/^Ao), and ga = 1, 1; 1, 2 to _Ro < l/(-ffo\/^Xvo)- In the 
latter case there are no solutions for a certain range of values of gi as described by Eq. 11411 . The thick 
lines show the (u = 0)-branch of R^/Ro, and the thin lines the {u — — l)-branch, respectively. 



For initial values Rq and mass parameters qi, that do not satisfy the above existence and stability 
conditions, the fate of the universe will be "big crunch"-like or "big rip"-like^, respectively. With these 
notations we mean that the scale factor a{t) or one of its derivatives H, q becomes singular in a finite 
time in the future. There is one exception, that may occur when qi and (72 are both positive, where the 
function F and thus the Hubble scale H approach constant values, while the horizon radius R goes to 
infinity, see Fig. 13 Another property of the cosmic evolution for negative qi is that no big crunch may 
occur. This can be seen from the time-derivative of the function KqF, 



which gets positive for H < (big crunch), thus preventing a further decrease of K^F and the final 
collapse of the scale factor. 



Up to this section, we analysed the evolution equation H1U|I for the scale factor analytically. Unfortunately, 
finding explicit solutions seems to be rather difficult because of the strongly non-linear form of the 
equation. Therefore, we solve it numerically, thereby realising that the form given by Eq. is not 
directly suitable for a numerical study because of the integral over the time t in the radius function R{t). 
This integral can be removed by differentiating with respect to t, leading to an ordinary differential 
equation of the order three, 




5 Numerical solutions 



a 



a 



+ 



(g2 - qi)KoF 
[qi - q2K^FY 



■ exp 



qi - q2KQF 



KoF-l 




^Recent analyses of future singularities and their properties can be found in Refs. [15]. 



which can be solved for the scale factor a{t) numerically. Note that one has to check whether the numerical 
solutions are also solutions of the original equation. This is not guaranteed since differentiating Eq. H1U|) 
possibly changes its set of solutions. Indeed, we encounter such "false" solutions in some cases. 

Regarding recent observations, we get acceptable solutions only for \qi\ to be of the order 1, which means 
that the relevant mass scale m should be near Aq^'' ^ 10^^ eV. Actually, the only known particles with 
such a low mass are neutrinos. This indicates that the influence of higher mass fields is suppressed, or 
these fields have decoupled, respectively. Unfortunately, the simple form of the RGEs Q and Q cannot 
account for a decoupling mechanism. Therefore, we assume in this work that the mass scale m is low 
enough today, so that the solutions are compatible with observations. However, note that at earlier 
cosmic times high-mass fields m should be relevant. 

Concerning the differential equation, we can fix several initial conditions by using observational results. 
These are today's value of the Hubble scale Hq = -{to) and the relative vacuum energy density Hao- 
Neglecting the spatial curvature {flko = — fc/d^ = 0) and considering only dust (with an equation of state 
parameter Q = 0, 5) and the CC as relevant energy forms in the present-day universe, the acceleration 
parameter go = f?(^o) is determined by Eq. q = + Q) + Q{Vlk — 1). Today's value of the horizon 
radius i?o is unknown, so we have to estimate it. Since it should be the largest physical length scale and 
the universe seems to be almost de Sitter-like, we assume the horizon radius to be around _Ro ~ Ij 2 • Hf^^ . 

For the numerical treatment every dimensionful quantity is expressed in terms of today's Hubble scale 
(Hubble units). In a ACDM universe with constant A and G, the cosmic age is denoted by to. In our 
calculations we used the following numbers from Ref. [16]: 

77o = 1,5 • lO^-^^QgV, r2Ao = 0,73, to 13, 7Gyr = 0, 99 • iJo"\ 
Ao = 2,98- lO^-^^GeV^, a| w 2, 34 • 10"^ eV, Go (1, 22 • 10^^ GeV)"^. 

The first observation from the numerical solutions is, that for a positive value of qi the cosmic age 
decreases with respect to the age to of the standard ACDM universe, whereas for a negative qi the age 
increases. Furthermore, the usually small value of 92 leads to a negligible time variation of Newton's 
constant G(t). 

To show the characteristic future cosmic evolution, we investigate four cases in more detail, which result 
from the parameter choices qi — ±2 and (72 — ±0, 1. Note that due to the suppression by the Planck 
scale, the realistic value of 52 should be much lower than ±0, 1. Here, we used a large value for q2 to show 
the differences due to the sign of (72 • Figures show the numerical results for different values of the 
initial radius i?o of the event horizon. The graphs in each of the four figures illustrate the scale factor a(t), 
the Hubble scale -ff(t), the acceleration q{t), the event horizon radius -R(t), and F{t) as functions of the 
cosmic time t, respectively. The last graph displays KqF as a function of the radius R/Rq. 

In section 0] we discussed the evolution equation analytically, and we found several conditions for the 
existence of stable de Sitter final states. These properties are also shown by the the numerical results. 
For negative values of qi, we observe only "big rip"- like solutions and de Sitter final states, whereas for 
positive qi, a big crunch may also occur, and all de Sitter states are unstable. Note that the catastrophic 
events, the "big rip" and the "big crunch", usually involve a high gravitational strength, implying that our 
calculations may not be reliable near these singular points. For positive values of qi and q2 (see Fig.CJ, 
we have not observed any "big rip"-like solutions. Then the final state may be either a "big crunch" or 
a forever expanding universe, where the Hubble scale approaches a finite positive value, but the event 
horizon radius R goes to infinity. This is a contradiction, because an asymptotically constant Hubble 
scale H > implies a finite event horizon radius R w in the far future, which is not the case here. 
Obviously, this numerical solution is not a solution of the original equation (|l()|l . For gi > and 52 < 
(see Fig. EJ the "big rip"-like events in the numerical solutions occur at a finite and large value of the 
horizon radius R. Again, this behaviour is not compatible with the vanishing of the horizon radius at 
such an event. Therefore, we can reject these numerical solutions, too. 

6 Conclusions 

We investigated a cosmological model, where the CC and the NC are determined by renormalisation 
group equations, which emerge from QFT on curved space-time. By choosing the inverse radius of the 
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cosmological event horizon as the renormalisation scale, we get time-dependent constants. Because of this, 
the evolution equation for the cosmic scale factor becomes more complicated than in standard ACDM- 
cosmology, leading to cosmological solutions with several very different future final states of the universe. 
We found "big crunch"-like and "big rip"-like solutions, but also stable de Sitter final states. Which of 
these states will be realised depends crucially on the field masses in the renormalisation group equations, 
and on the initial value of today's radius of the event horizon. In this context, we derived some conditions 
on the existence of stable dc Sitter states. Furthermore, the cosmic evolution was analysed numerically 
for different values of the field masses. For a realistic cosmic behaviour, we have to require that the masses 
in the RGEs should be quite low, implying the need for some suppression or decoupling mechanism for 
high-mass fields. However, such a mechanism for the CC and the NC has not been found yet [12]. This 
also restricts the main focus of this work to the future behaviour of the universe, because for the cosmic 
evolution at early times quantum fields with high masses should be taken into account. Moreover, the 
regimes with a high gravitational strength need a deeper investigation, because such conditions arc given 
not only at early times, but also at the singularities in the future. Finally, wc> ('onc;lude, that the specific 
choice of the renormalisation scale in this work leads to cosmological solutions, that may become singular 
in a finite time without introducing exotic forms of matter. 
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Figure 4: The cosmological evolution for the parameter choice qi — —2 and 52 = —0, 1 and different values of 
the initial event horizon radius _Ro- The fate of this universe is either a stable de Sitter state when 
choosing Ro = 1, 20; 1, 30, or a big rip (BR) in the case _Ro = 1, 10; 1, 15. K^F is bounded from above. 
Nomenclature: Scale factor a, Hubble scale H — acceleration q — event horizon radius R and its 
initial value Ro. For the function KqF see Eqs. (1121, 0, and jnj, for the mass parameters qi , g2 see 
Eqs. ^ and 
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Figure 5: For the parameter choice q\ = —2 and = +0, 1 the cosmic evolution is not very different from the 
case §2 = —0, 1 (Fig. [IJ. In the future, there is either a stable de Sitter state for _Ro = 1, 20; 1, 30, or a 
big rip (BR) when iio = 1, 10; 1, 15. KqF is bounded from below. Nomenclature; Scale factor a, Hubble 
scale H = ^, acceleration q = event horizon radius R and its initial value Rq. For the function KqF 
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